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ABSTRACT 


Certain  jobs  in  large  projects  do  not  have  a  unique 
"location"  in  the  critical  path  network;  they  may  be  moved 
into  certain  slack  intervals,  for  example,  or  may  even  be 
divisible  into  smaller  subtasks,  and  "tucked  in"  at  several 
locations.  An  example  is  the  painting  activity  during  air¬ 
craft  maintenance  and  overhaul  which  may  logically  be 
performed  at  several  points  during  the  overhaul  cycle. 

In  a  previous  paper  £4],  an  analysis  was  given  of  a  model  in 
which  a  single  job  can  be  divided  up  in  any  manner  among  an 
arbitrary  number  of  locations;  the  resulting  algorithm  was 
of  the  optimal  network  flow  type,  which  can  be  simply  and 
efficiently  solved  using  available  computer  codes.  In  the 
first  part  of  the  present  paper,  this  model  is  extended  to 
multiple  jobs  of  divisible  type.  The  general  approach  is 
via  the  decomposition  method  of  linear  programming;  however, 
the  resulting  algorithm  is  again  fairly  simple.  Optimal 
cost-time  solutions,  possibly  infinite,  (or  optimal  network 
flow  solutions,  possibly  infeasible)  are  generated  by  known 
algorithms.  The  resulting  schedules  or  cuts  are  then 
combined  in  a  simple,  special-structure  linear  program  whose 
dimensionality  is  equal  to  the  number  of  divisible  groups. 
Bounds  on  the  nonintegrality  of  the  final  allocations  can 
also  be  determined. 

When  these  special  jobs  can  only  be  moved  about  the  network 
in  their  entirety,  or  in  certain  indivisible  modules,  the 
problem  takes  on  the  form  of  an  integer  program.  In  the 
second  part  of  the  paper,  a  branch-and--bound  procedure  will 
be  given  for  the  problem  of  movable  activities,  together  with 
efficient  heuristics  for  arbitrating  and  bounding  these 
locations,  using  only  the  ordinary  critical-path  algorithm. 

Examples  are  given  for  both  models. 


DIVISIBLE  AND  MOVABLE  ACTIVITIES  IN  CRITICAL  PATH  ANALYSIS 


by 

William  S.  Jewell 


0.  INTRODUCTION 

Certain  jobs  In  large  projects  do  not  have  a  unique  "location"  In  the 
critical  path  network;  they  may  be  moved  into  certain  slack  intervals,  for  example, 
or  may  even  be  divisible  into  smaller  subtasks,  and  "tucked  in"  at  several 
locations.  An  example  is  the  painting  activity  during  aircraft  maintenance  and 
overhaul  which  may  logically  be  performed  at  several  points  during  the  overhaul 
cycle. 

In  a  previous  paper  [4],  an  analysis  was  given  of  a  model  in  which  a  single 
job  can  be  divided  up  in  any  manner  among  an  arbitrary  number  of  locations;  the 
resulting  algorithm  was  of  the  optimal  network  flow  type,  which  can  be  simply  and 
efficiently  solved  using  available  computer  codes.  In  the  first  part  of  the 
present  paper,  this  model  is  extended  to  multiple  Jobs  of  divisible  type  giving  a 
decomposition  algorithm  with  special  structure  in  both  the  master  and  subprograms. 

Another  possibility  is  that  such  a  task  cannot  be  subdivided  as  finely  as 
desired,  but  can  only  be  moved  around  in  certain  modules,  or  in  its  entirety.  In 
the  second  part  of  this  paper,  a  branch-and-bound  algorithm  is  developed  which 
uses  the  ordinary  critical  path  algorithm  in  an  efficient  manner. 
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1.  FORMULATION 


Assume  that  a  project  network  of  A  arcs  ( activities )  and  N  nodes 
(events)  is  given.  Arc  ij  represents  activity  beginning  at  node  i  and 
ending  at  node  j  ,  nodes  1  and  N  represent  Start  and  Finish  of  the  project, 
respectively.  Activity  durations  {T^}  are  given  for  each  arc,  and  the  problem 
is  to  select  potentials  (event  epochs)  {v^}  for  each  node  which  minimize  the 
total  project  duration,  F  ”  vn  ~  vl  *  As  usual,  we  assume  all  finite, 

that  no  directed  cycle  exists  with  positive  sum  of  ,  and  at  least  one 

directed  path  exists  from  node  1  to  N  . 

Additionally,  we  assume  there  are  D  divisible  activities  which  require 
duration  ,  (d  «  1,2,  . . . ,  D)  .  These  durations  may  be  divided  up  in 
arbitrary  amounts  and  carried  out  at  any  or  all  of  a  certain  number  of  "locations" 
in  the  network,  specified  in  advance  for  each  such  activity.  For  simplicity  in 
notation,  artificial  arcs  are  added,  if  necessary,  so  that  an  arc  is  either 
a  location  for  a  single  divisible  task,  or  for  none.  In  other  words,  partition 

the  class  of  all  arcs,  A  ,  into  D  +  1  exclusive  subsets,  Aq  ,  A^ . Ap  , 

and  read  "ij  e  A^"  as  "activity  ij  belongs  to  divisible  class  d";  conversely, 
"d  e  A  ^(ij)"  reads  "d  is  the  (unique)  class  to  which  ij  belongs."  Let  A. 


be  the  number  of  arcs  in  class  d  ,  so  that  A^  +  A^  +  • • .  +  A^  «  A  . 

Finally,  define  t^j  as  the  additional  duration  expended  on  divisible  work 
at  location  ij  ;  the  initial  primal  problem  is  then: 


(la) 

Min  F  “  vjj  ”  vi 

(lb) 

vj  ‘  vi  ‘  fcij 

■n 

H 

A  1 

V  ij  e  A 

(lc) 

‘ij 

>  0 

(Id) 

l  tu 

ijeA.  « 

-Td 

d  ■  1,2,  ....  D 

3 


It  is  immaterial  whether  we  set  T  «  0  ,  or  drop  the  t  .  from  arcs  in  A 

o  ij  o 

The  constraint  matrix  for  the  example  of  Section  8  is  shown  in  Figure  1.  The  upper 
part  of  the  matrix  is  the  same  as  that  of  a  cost-time  critical  path  (CPM)  problem 
[3],  [4]  where  the  objective  is,  typically: 


(la’) 


Min  F 


Q(v„  "  v  )  -  l  C  t  +  fixed  terms 
N  1  ij  eA  iJ  lj 


Here  is  a  unit  savings  gotten  by  "lengthening"  ij  by  an  amount  t^j  above 

the  "crash"  duration  T^  ;  usually  there  is  an  upper  bound  on  t^  to  limit  paid 
lengthening  to  some  "normal"  duration.  Q  is  a  Lagrange  multiplier  varied  during 
the  course  of  the  cost-time  algorithm,  which  has  a  dual  interpretation  of 
exogenous  flow. 

The  constraints  (Id),  on  the  other  hand,  are  like  the  "bundle  constraints" 
of  multi-commodity  flow  problems  [3],  (6),  [9],  except  that  here  they  are  adjoined 
to  the  transpose  of  the  usual  Kirahoff  flow  matrix. 

The  dual  initial  problem  is: 


(2a) 

(2b) 

(2c) 

(2d) 


Max 


l-  I 

ij  eA 


Vu 


D 

l 

d-1 


Vd 


Xij) 


i  -  1 
i  -  N 
otherwise 


xij  2  0 


‘n  -  y  -1 

iJ  A  1Cij) 


ij  e  A 


which  is  the  usual  dual  longest-zoute  problem,  except  for  a  possible  profit  for 
simutaneously  increasing  several  lower  bounds,  rather  like  an  optimal  capacity 
contraction  problem. 

It  is  the  confluence  of  these  different,  but  special  models  that  r  -kes  the 
analysis  of  (1)  interesting  from  a  pedagogical  point  of  view,  providing  we  can 
retain  the  computational  simplicity  of  these  special  structures. 
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2.  DECOMPOSITION 

We  rhall  solve  (1)  by  using  the  decomposition  method  ( see  [2]  or  [10, 

Chapter  10]),  with  (Id)  as  the  master  problem  and  (lb),  (lc)  as  the  subproblem. 

The  master  will  be  solved  by  the  simplex  method  (with  some  simplifications),  and 
the  subproblera  will  be  solved  by  special  CPM  or  optimal  flow  algorithms.  For 
convenience,  we  review  the  principles  of  this  method  in  our  notation. 

The  set  of  all  solutions  to  (lb) ,  (lc)  is  a  convex  polytope  with  both 
finite  and  unbounded  regions;  in  fact,  the  unbounded  region  usually  turns  out 

to  be  important  in  the  application  of  the  algorithm.  Let  (k  "  1>2 . K) 

be  the  set  of  all  extremal  solutions  to  the  finite  region  (a  polyhedron)  and 
{wi'uij}  ^  c  • • • *  L)  be  the  extremal  rayB  bounding  the  infinite  part  of  the 

solution  space  (a  convex  polyhedral  cone) .  From  a  well-known  theorem  of  convex 
sets  [10],  any  feasible  set  of  solutions  to  (lb),  (lc)  can  be  represented 

as  follows: 
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(5a)  Min  F  =  £  p\  +  £  gV 
k  K  l 


(5b) 

k  * 

-  1 

(5c) 

l  Td\  +  l  Ud9i 

k  l 

-  Td  (d  -  1,2, 

....  D) 

(5d) 

\  >  0  (k  -  1,2, 

...,  K)  ;  8  >  0  (i  -  1,2, 

....  L) 

Thus,  the  problem  is  changed  to  one  with  D  +  1  constraints  but  many  variables. 
The  new  form  of  the  dual  master  is: 


(6a) 

Max  L  »  a  +  2 
d 

(6b) 

0  +  t  Td*d 

d 

V  1 

(k  - 

1,2, 

....  K) 

(6c) 

d 

<G* 

a  - 

1,2, 

....  L) 

(6d) 

a  ,  *d 

unrestricted 

(d  - 

1,2, 

.  ...  D) 

o  and  the  will  always  be  nonnegative  in  any  case  of  interest,  since  weakenip. 

(5b),  (5c)  to  inequalities  will  not  change  any  optimal  solution. 

As  is  usual  in  decomposition,  one  does  not  need  all  extremal  solutions  and 
rays,  to  start  with,  but  only  the  D  +  1  necessary  to  form  an  initial  feasible 
basis  to  (5).  New  candidates  will  be  generated  as  needed  by  solving  the  subproblem 
and  added  during  the  course  of  the  algorithm.  At  any  stage,  if  the  current  optimal 
dual  variables  are  »  a  finite  candidate  solution  k  «  *  to  the  subproblem 

can  be  checked  by  testing  dual  constraint  (6b).  Thus,  if: 


A  r  o  * 

F  -  l  Vd 


<  o 


(7) 


6 


(7') 


7 

then  this  optimal  extremal  solution  c^ear^y  a  candidate  for  the 

master.  But  if  the  minimal  solution  does  not  satisfy  (7*),  then  no  new  candidates 
exist,  and  the  current  master  must  be  optimal.  Note  that  any  unbounded  optimal 
solution  to  (8)  will  always  satisfy  (7*)« 
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3.  THE  NATURE  OF  UNBOUNDED  SOLUTIONS 


The  associated  subproblem  (8),  (9)  is  in  cost-time  problem  form,  except 
there  is  no  limit  on  the  savings  possible  for  extending  activities  with  large 
C^j  ,  i.e.,  the  optimal  value  of  E  may  be  unbounded.  In  dual  terms,  the 
exogenous  flow  Q  -  1  may  be  insufficient  to  "cover"  several  large  lower  bounds 


on  flow  in  arcs  which  are  in  a  parallel  configuration,  i.e.,  (9)  may  be 
infeasible . 


Figure  2  clarifies  the  nature  of  unbounded  solutions.  Suppose  the  set  of 
all  nodes,  N  ,  is  partitioned  into  a  subset  X  containing  node  1  ,  and  the 
complementary  set  N  -  X  containing  node  N  ,  the  set  of  all  arcs  with  one 
extremity  in  X  ,  the  other  in  W  -  X  ,  is  called  a  cut.  In  Figure  2,  we  have 
an  oriented  cut  X  ,  which  all  arcs  ij  are  in  (X,N  -  X)  ,  i.e.,  they  pass 
"from  left  to  right."  Most  of  the  cuts  in  actual  critical  path  networks  are  so 


oriented. 


{v'ij} 


be  an  (8)  feasible  solution  to  Figure  2.  Then  it  is  clear 


(ID 


v^  ■  v°  +  0  ieN-X 

‘ij  ’  ‘ij  +  9  ij  e  (X’W  ”  X) 
t|j  *  t°j  otherwise 


is  also  feasible  for  all  0  <  0  <  “  .  Thus,  an  oriented  cut  represents  a 
possible  unbounded  solution  to  (8).  One  can  easily  show  that  no  other  possible 
infinite  rays  exists,  assuming  that  the  usual  ambiguity  in  the  {v^}  is  resolved 
by,  say,  setting  v^  -  0  .  For  example,  reversing  a  single  arc  in  Figure  2 
leads  immediately  to  an  upper  bound  on  6  . 


V1  V2  V3  V4  V5  C12  C13  fc14  C23  *24  C35  C45 


FIGURE  1: 


CONSTRAINT  MATRIX  AND  REQUIREMENTS  VECTOR  FOR  EXAMPLE 
OF  SECTION  8 


FIGURE  2:  ORIENTED  CUT  OBTAINED  FOR  UNBOUNDED  SOLUTIONS  TO  SUBPROBLEM 
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Let  a.(X)  be  the  number  of  arcs  of  divisible  class  d  across  oriented  cut  X  . 
d 

Unbounded  solutions  to  (8)  can  be  found  by  looking  for  infeasible  solutions  to 
(9b),  (9c).  For  example,  we  could  set  the  exogenous  flow  equal  to  Q  ,  and  solve 
a  minimum  flow  problem;  if  Q  »  Min  Q  >  1  ,  then  (9b),  (9c)  is  infeasible. 
However,  by  using  the  min  flow  =  max  out  theorem,  it  follows  that  an  oriented 
cut  X  represents  an  unbounded  solution  to  (8)  iff: 

(12)  l  ad(X)rr°  >  1  . 

d 

An  economic  interpretation  of  (12)  is  that  unbounded  solutions  arise  when  the 
total  unit  profitability  for  expanding  an  oriented  cut  is  greater  than  the 
resulting  "cost"  of  increasing  the  project  duration  by  one  time  unit. 
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4.  THE  MASTER  PROBLEM 

The  master  problem  needs  at  least  D  +  1  independent  solutions  to  begin. 

For  example,  if  an  extremal  solution  k  to  (lb),  (lc)  (or  even  just  a  feasible 
solution)  is  known  then  the  associated  column  in  (S)  is,  transposed: 

<l3>  [*«S4 . TS]  • 

(the  first  position  is  for  the  objective  function).  An  efficient  bounded  starting 
solution  is  given  in  Section  5. 

However,  in  many  problems,  the  rays  associated  with  unbounded  solutions  play 
an  important  role.  From  the  above  discussion,  the  candidate  column  associated  with 
an  oriented  cut  X  will  be,  transposed: 

(14)  [l;0;a1(X),a2(X),  ....  aQ(X)]  . 

Since  these  oriented  cuts  and  their  associated  coefficients  can  be  easily 
generated  by  various  simple  algorithms,  it  is  advantageous  to  use  as  many  of  them 
as  necessary  in  the  initial  basis;  indeed,  the  optimal  basis  often  consists  of 
D  rays  and  one  finite  solution.  For  small  problems,  one  can  even  enumerate  all 
oriented  cuts  and  start  with  the  "best"  ones  (largest  sums  7  a. (X)). 

The  resulting  master  program  is  solved  by  the  simplex  method,  and  then 
resolved  as  new  columns  are  generated  by  the  subproblem  until  finally  a  finite 
solution  is  generated  which  fails  test  (7'),  It  is  advantageous  to  carry  along 
an  adjoined  unit  matrix,  so  that  the  current  inverse  of  the  basis  can  be  used  to 
update  encering  columns.  Columns  dropped  from  the  basis  tend  not  to  reappear, 
but  there  is  no  guarantee  of  this. 


13 


5.  AN  EFFICIENT  BOUNDED  STARTING  SOLUTION 

It  is  inefficient  to  use  the  ordinary  critical  path  solution  to  (1) 

(with  all  T,  *  0)  as  a  starting  solution,  since  even  if  v„  -  v.  is  maintained 
a  NX 

at  its  minimal  duration,  F°  ,  there  is  almost  always  enough  slack  (float) 
off  the  critical  path  to  "tuck  in"  some  divisible  time  "for  free."  The  following 
procedure  is  suggested  as  a  good  solution  to  use  as  the  k  »  1  column  in  the  master. 
(It  is  assumed  that  the  critical  path  problem  can  be  solved  repeatedly  in  an 
efficient  manner,  and  that  a  ranked  list  of  currently  unallocated  T^  is  kept, 
and  continuously  updated.) 


(15)  (a)  Set  all  T,  *■  0  ,  and  solve  (1)  for  its  minimal  solution 
a 

|F°jv°;t°j  «  oj  ,  with  earliest  possible  event  times.  The  resulting 
tight  activities  constitute  an  early  tree  with  a  trunk  (critical  path) 
from  node  1  to  N  ,  and  branches  to  all  other  nodes .  For  all  nodes 
on  the  trunk,  v*  «  v°  ,  (F1  «*  F°)  ,  and  for  all  arcs  on  the  critical 
path  t^  «  0  . 


(b)  If  ij  is  a  aotree  arc  belonging  to  a  currently  unallocated  class 

A  ^(ij)  *  d  ,  let  t°  <=  Min  fv°  -  v°  -  T  -  t°  ;  currently  unallocated 
ij  L  J  1  1J 

tJ  .  (Arcs  may  be  taken  in  any  order.) 

(c)  Keeping  the  t°j  of  (b)  fixed,  resolve  for  the  late  tree  (latest  possible 
event  times).  Repeat  (b)  for  any  cotree  arcs  in  a  currently  unallocated 
cohort . 

(d)  Keeping  the  current  t°^  fixed,  find  the  earliest  |v^|  and  latest 
{Vi}  event  times  from  two  passes  of  the  critical  path  algorithm.  For 
each  currently  unallocated  class  d  ,  rank  the  arcs  ij  c  A^  on  the 
basis  of  largest  available  slack 


8ij 
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(e)  Select  the  class  d  with  current!/  largest  unallocated  and 

positive  slack  in  some  arc.  For  the  ij  e  with  largest  , 

let  t°j  <=  Min  [s^j  ;  currently  unallocated  T^] . 

(f)  Repeat  (d)  finding  the  new  event  times  until  either: 

(i)  All  {T,}  are  allocated;  the  solution  is  trivial. 

d 

or, 

(ii)  There  exist  unallocated  divisible  classes,  and  all  arcs 

in  these  classes  have  zero  slack,  i.e.,  there  are  multiple 
critical  paths  of  tight  arcs.  Set  t^  >=  t°^  ,  and  the 
v^  to  either  of  the  values  obtained  during  the  last  pass 
(d)  -  (f)  during  which  no  allocations  occurred. 


There  are,  of  course,  many  other  good  heuristic  procedures  which  can  be  used. 
The  advantage  of  the  above  method  is  that  it  requires  only  two  critical  path 
solutions  for  (a)  and  (b),  which  often  takes  care  of  much  of  the  "tuckable" 


t^  ,  and  then  only  two  solutions  plus  some  ranking  per  pass  (d)  -  (f).  Since  (e) 
either  increases  the  number  of  tight  arcs,  or  completes  the  allocation  of  some 
class,  the  procedure  is  finite. 

We  believe  this  process  gives  a  reasonably  good  starting  solution.  However, 
notice  that  it  need  not  be  near  the  final  optimal  allocation  (which  may  require 
large  allocations  on  one  arc),  nor  is  it  optimal  in  the  sense  of  allocating  the 


maximal  possible  sum  of  all 


since  fractional  solutions  are  not  allowed. 


Once  the  partial  network  of  tight  arcs  obtained  in  this  solution  has  been 


identified,  it  is  trivial  to  devise  algorithms  to  generate  the  D  oriented  cuts 


needed  for  the  rest  of  the  initial  basis. 
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6.  SOLVING  THE  SUBPROBLEM 

Some  remarks  on  solving  (8),  (9)  are  in  order.  If  a  cost-time  (CPM) 
computer  code  [3],  [7]  is  available,  then  (8)  could  be  solved  with  functional 


(8a') 


Min  E  =  Q(v 


N 


1  ijeA 


and 


(15)  (crash  duration)^  =  T^,  , 

(16)  (normal  duration)^  “  7^  +  » 

where  is  some  large  but  finite  number,  say  »  M  for  all  arcs.  To 

find  infinite  rays,  it  would  be  best  if  M  were  set  larger  than  any  number  to 
which  it  were  compared  during  the  course  of  calculation;  and  it  was  possible  to 
run  the  algorithm  backwards,  i.e.,  with  decreasing  Q  .  An  infeasible  solution  would 
be  detected  at  the  breakpoint  where  Q  attained  unity,  and  the  cut  would  be  the 
set  of  arcs  with  t^  «  M  . 

An  alternative  method  which  uses  the  CPM  code  in  the  usual  manner  Is  to  set: 


(17) 


Kii  “  T  -1 
ij  A  X(ij) 


V  ij  A  -  A  . 

o 


This  will  give  only  finite  solutions  to  the  master  problem.  If  there  are  an 
infinite  number  of  optimal  solutions  to  Q  “  1  (i.e.,  there  are  two  breakpoints 
Q  *  a  and  Q  ■  8  ,  with  a  <  1  <  6),  then  both  extreme  solutions  should  be 
furnished  to  the  master  for  optimal  mixing. 

Another  possibility  is  to  work  with  the  dual  (9),  using  an  out-of-kilter 
code  [3J.  Figure  3  shows  the  complementary  slackness  diagram  which  is  appropriate 
to  the  dui]  flow  problem,  assuming  the  code  is  of  minimizing  type.  All  other 
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things  being  equal,  it  is  desirable  to  have  t^j  as  large  as  possible,  so  that  a 
feasible  solution  x^  »  should  be  interpreted  as  being  a  corner  solution 

(i.e.,  the  relation  (8b)  is  tight).  If  the  dual  is  infeasible,  the  labelling 
subroutine  of  the  out-of-kilter  method  will  stop  with  some  nodes  labelled,  and 
the  remainder  not;  these  labels  then  define  the  desired  oriented  cut  X  . 

Alternatively,  one  can  place  two  arcs  in  parallel,  with  zero  lower  flow 
bounds,  and  use  parameters: 


(18) 


Arc  (ij)' 
Arc  (ij)" 


giving  the  composite  complementary  slackness  diagram  shown  by  the  "staircase" 
dotted  lines  in  Figure  3.  This  procedure  gives  only  finite  solutions  to  the 
master.  Again,  if  solutions  on  the  "riser"  of  Figure  3  are  obtained,  both 
endpoints  should  be  given  to  the  master. 


Most  computer  codes  require  integer  ,  whereas  ours  may  be  fractional. 

However,  there  alvzays  exists  some  large  integer  J  such  that  “  JC^ 

integer  for  all  ij  .  Hence,  an  equivalent  problem  can  be  run  with  integer  , 

and  flow  and  breakpoint  0  ■  J  .  The  resulting  flows  and  functionals  will  all 


be  too  large  by  the  same  factor;  however,  the  new  plan  or  cut  will  still  be  the  same. 


Section  9  discusses  a  choice  of  J  which  need  not  be  changed  at  every  pass  of 


the  subproblem. 


_______  Infeasible  Solutions  rc:-?*ble 

____  Feasible  Solutions  Only 

FIGURE  3:  COMPLEMENTARY  SLACKNESS  DIAGRAM  FOR  SOLVING  SUBPROBLEM 
WITH  MINIMAL  COST  OUT-OF-KILTER  ALGORITHM 


18 


7 .  THE  0 IVISLBLE  ACTIVITIES  ALGORITHM 

Figure  4  summarizes,  in  flow  chart  form,  the  final  decomposition  algorithm 
proposed  for  solving  the  divisible  activities  problem.  Here  t  is  the  index 
of  iterations,  t  -  0,1,2,  ...  Notice  that  the  problem  is  always  feasible 
(assuming  the  ordinary  critical  path  problem  is)  and  bounded. 


FIGURE  4:  FLOW  DIAGRAM  OF  DIVISIBLE  ACTIVITIES  ALGORITHM 
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8.  EXAMPLES  OF  DIVISIBLE  ACTIVITY  ALLOCATION 

Figure  5  shows  the  network  corresponding  to  the  matrix  of  Figure  1. 

The  numbers  on  the  nodes  are  the  node  indices;  the  numbers  on  the  arcs  represent 

the  }  .  Figure  6  shows  the  early  tree  solution  when  *  0  and 

corresponding  event  times  on  the  nodes,  the  late  solution  also  has  F  »  8  , 

but  v*1  »  5  . 

4 

Henceforth,  the  node  numbers  will  represent  the  {v^}  »  anc*  the  arc  numbers 
will  be  the  {t^J  .  Tight  arcs  will  be  solid  and  slack  arcs  dotted.  Figure  7 
shows  the  initial  b  unded  solution  obtained  when 


Ti;Ti 


3  ;  t2  >  T* 


2 


and  the  algorithm  in  (15)  is  used.  All  arcs  are  tight. 


1.  For  the  first  example,  suppose  T^  *  4  and  T^  *  3  ,  and  choose  as  cuts 

X1  =  {1}  ;  X2  =  (1,2, 3, 4}  . 

The  initial  form  of  the  master  tableau  is: 


a 

"i 
*2 

where  the  top  row  is  the  functional  F  ,  and  a  unit  matrix  has  been  adjoined  to 
provide  the  basis  inverse. 

The  usual  reduction  methods  give  the  reduced  tableau: 


8 

B 

B 

0 

0 

0 

0 

1 

0 

D 

1 

0 

0 

3 

3 

0 

H 

0 

1 

0 

2 

2 

3 

0 

0 

1 
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Ti  »  ij  -  3 

T2  -  *2  '  2 
Arc  numbers  - 

FOR  FIRST  EXAMPLE, 


FIGURE  7:  INITIAL  BOUNDED  SOLUTION  (1  ) 

WITH  UNALLOCATED  T,  T  1 
1  ’  l2 


,5)} 
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from  which  we  conclude  X^  »  1  ,  0^  »  1/3  ,  ©2  *  1/2  ;  o  *  6  ,  ir^  *  1/3  ,  =  1/2 

and  F  =  8  5/6  .  The  resulting  allocation  after  the  first  iteration  is  shown 
in  Figure  8. 

Next  a  flow  subproblem  with  C  *  1/3  (ij  e  A^)  and  C  =  1/2  (ij  e  A^) 

is  attempted,  and  an  infeasible  solution  obtained.  A  flow  Q  “  1  2/3  will 

"cov-'-r"  all  the  arcs,  but  attempts  to  reduce  Q  to  1  lead  to  a  max  cut 

3 

X  =  {1,2}  .  Consequently,  a  column  3^  is  added  to  the  master  with  values 
[1,0, 2, 2]  transPose  f  after  updating  in  the  usual  way.  0^  drops  from  the  basis 
and  the  new  optimal  solution  is  X^  «=  1  ,  0^  c  0+  ,  6^  =  1/2  ;  o  =  7  ,  *  0  , 

if 2  “  1/2  and  F  =  8  1/2  .  The  second  allocation  is  shown  in  Figure  9. 

When  the  flow  subproblem  is  solved,  a  finite  flow  solution  is  obtained,  with 
X12  =  x23  =  ^  ’  X34  =  X35  “  x45  =  1/2  and  K  =  7  1/2  .  The  corresponding  primal 

would  be  t^j  ~  0  t  ^3  15  1  *  ^14  ”  ^  *  ^23  *  ^  *  ^24  =  ^  *  ^35  *  ®  *  ^45  =  1  ^  ® 

(0  <  0  <  <»)  ,  but  we  know  already  from  duality  that  E  =  7  1/2  (V  0)  ,  test  (71) 
fails,  and  this  is  not  a  candidate  solution.  Hence,  Figure  9  is  the  optimal 
allocation  to  the  first  example. 

II.  As  a  second  example,  take  «  5  and  T2  «  3  .  The  same  sequence  of 

tableaux  is  obtained  except  that  0^  displaces  02  ,  and  X^  »  1  ,  0^  «  1/3  , 

9^  “  1/2  ;  a  »  6  2/3  ,  «  1/2  ,  tt ^  “  6  (Figure  10).  The  optimal  flow  dual 

subproblem  has  x^  *  x^  -  1/2  and  K  «  6  1/3  .  Since  test  (7')  is  satisfied, 


cuts  currently 
in  use 


FIGURE  8:  INITIAL  MASTER  SOLUTION  FOR  FIRST  EXAMPLE,  WITH  ALL 
DIVISIBLE  TIME  ALLOCATED 


FIGURE  9:  SECOND  (AND  OPTIMAL)  SOLUTION  TO  FIRST 
EXAMPLE 


FIGURE  10:  SECOND  MASTER  SOLUTION  TO  SECOND  EXAMFLE 
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we  must  find  the  new  finite  primal  subproblem  (X^)  ,  which  is  shown  in  Figure  11. 
Note  that  fractional  extremal  solutions  to  subproblems  can  arise  (with  fractional 
,  and  that  (3,4)  is  slack  for  the  first  time.  A  new  column 
[8  l/2;l;5,3]trans^ose,  is  added  to  the  master  after  updating,  and  X£  drives  out 
(say)  X^  ,  leaving  8^  ■  8^  at  zero  level  in  the  basis,  o  ■  6  1/3  ,  and 
ir^  *=  1/3  ,  ^2  =  1/6  as  before,  and  this  time  the  test  fails  as  an  equality.  In 
other  words,  the  bounded  extremal  solution  X2  “  1  is  by  itself  the  optimal 
solution  (another  example  of  this  occurs  in  the  example  given  in  [4],  where 
D  *  1)  .  The  additional  unit  of  T^  is  handled  without  increasing  the  project 
duration. 

III.  As  the  third  example,  let  T^  *  4  and  T2  “  2  .  By  inspection  of  X^  we 
see  that  all  is  taken  care  of,  and  the  cheapest  way  to  handle  the  remaining 

unit  of  is  to  put  as  much  of  it  in  parallel  as  possible,  i.e.,  to  divide  it 

into  thirds.  The  optimal  solution  is  shown  in  Figure  12,  where  X^  <*  1  , 

0j_  =  1/3  ,  83  -  0+  ;  a  =  6  2/3  ,  -  1/3  ,  tt2  - 


1/6  . 


0 


FIGURE  11:  NEW  FINITE  EXTREMAL  SOLUTION  U2> 
FOR  SECOND  EXAMPLE  (ALSO  OPTIMAL) 


FIGURE  12: 


OPTIMAL  SOLUTION  TO  THIRD  EXAMPLE 
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9.  SOLUTION  MODULARITY 

The  reason  that  fractional  solutions  arise  is  because,  when  parallel 

critical  arcs  are  available  for  allocation,  it  is  optimal  to  divide  each  hour  of 
unallocated  time  into  the  fractions  1/a^  »  and  spread  the  work  uniformly  over  the 
critical  paths.  Following  the  arguments  advanced  in  [3],  we  see  that  if  D  =  1  , 
then  only  solutions  in  multiples  of  1/1, 1/2,  ...,  1/A^  could  be  obtained.  With 
D  >  1  ,  it  follows  that,  if 

(19)  A*  -  max  A.  , 

d=l ,2 , . . .  ,D 

then  multiplication  of  the  original  duration  data  by  (A^) !  must  lead  to  a 
solution  in  integers.  In  the  terminology  of  [6],  the  problem  is  A^  -modular. 

This  result  is  primarily  of  theoretical  interest;  however,  it  does  indicate 
that  J  =  (A*)'  is  a  suitable  choice  when  using  an  integer  code  as  described  in 


Section  6. 


27 


10.  MOVABLE  ACTIVITIES 

A  different  type  of  problem  arises  when  the  constraints  (lc)  (Id)  are 
replaced  by 


(20) 


«Ld '« 

' 0  "  V<«,  v  1J  • 


In  this  case,  we  consider  that  the  work  of  type  d  is  not  divisible,  but  is 
only  movable ,  i.e.,  it  may  be  performed  only  at  one  of  several  possible  locations 
in  the  project  network.  This  problem  is  one  of  binary  programing,  since  one  can 
introduce  binary  variables,  d  ,  with 


(21) 


t.,  *  i./T  , 

ij  lj  A_1(ij) 


Similar  formulations  arise  if  we  consider  that 
indivisible  modules  of  work 


k’  |Tdl,Td2 . TdMd 


the  job  content  consists  of 
such  that 


“d 

l 

m=l 


dm 


(d  «  1,2,  ....  D)  . 


Conceptually,  one  separates  arcs  in  group  d  into  arcs  in  series,  and  treats 

each  module  as  a  separate  movable  group.  A  special  example  of  interest  is  when 
these  modules  are  of  the  same  length,  say  one  unit  of  time;  this  would  be  the 
problem  of  divisible  activities  iaith  integer  solutions . 

There  are  also  mixed-integer  formulations  possible,  but  these  will  be  seen  to 
be  a  special  case  of  the  following  approach.  They  are  also  of  less  practical 


interest . 
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11.  SEPARATION  IN  A  BRANCH-AND-BOUND  METHOD 

The  method  proposed  to  solve  the  problem  of  movable  activities  la  a  branch- 
and-bound  procedure,  modified  to  take  account  of  the  special  structure  of  the 
problem.  Readers  not  familiar  with  this  approach  are  referred  to  Lawler  and 
Wood  [8],  or  Bertier  and  Roy  [1]. 

The  first  element  to  be  considered  is  the  branching  (arbitration) ,  i.e.,  the 
separation  of  the  solution  space  Into  subspaces  where  certain  of  the  binary 
variables  are  fixed  once  and  for  all.  This  separation  is  conveniently  represented 
by  an  arbors scence ,  where  each  node  represents  a  possible  partition  of  the  solution 
space,  and  each  arc  shows  the  branching  between  a  given  arbitration  and  the  further 
possible  partitions. 

The  natural  integer  programming  separation  would  come  about  by  considering 

some  (possibly  adaptive)  ordering  (ijj^ji^Jj*^! 3*  •  ••)  the  6^  »  and  making 

the  arbitration  on  whether  6  =0  or  1  .  This  leads  to  the  arborescence  shown 

in  Figure  13.  However,  this  method  tends  to  postpone  the  labor  involved  in  making 

and  evaluating  the  decision  as  to  where  the  work  T,  should  be  performed,  since 

a 

the  bound  used  for  6^  =  0  always  leads  to  further  exploration  of  that  "branch," 
until  some  6  *  1  is  selected. 

A  more  compact  representation  is  the  task-oriented  arbitration  shown  in 

Figure  14,  in  which  a  (possibly  adaptive)  ordering  of  the  movable  tasks  is  made, 

and  each  separation  represents  a  definite  allocation  of  T,  to  one  of  A, 

a  d 

locations. 

We  call  an  arbitration  of  order  k  a  problem  in  which  exactly  k  cf  the 
have  been  fixed  at  0  or  1  .  The  result  can  be  considered  a  new  problem 


T'  - 

U 


T  +  T 

ij  A-X(ij> 


ij  arbitrated  and  6 


otherwise; 


No  decisions  made 


decided 


FIGURE  13:  ARC-ORIENTED  SEPARATION  OF  MOVABLE  ACTIVITY  SOLUTION  SPACE 


No  decisions  made 


FIGURE  14:  TASK-ORIENTED  SEPARATION  OF  MOVABLE  ACTIVITY  SOLUTION  SPACE 
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the  t^j  are  then  eliminated  from  the  arbitrated  arcs,  giving  a  reduced  problem 

in  A  -  A  -  k  free  t . .  . 
o  ij 

If  we  use  the  task  oriented  separation  of  Figure  4,  we  see  that  the  arbitration 
jumps  from  order  0  to  order  A,  from  the  first  to  the  second  level,  leav  .ng  a 


new  problem  with  one  group  of  t^  eliminated,  and  only  D  -  1  more  levels  to  go. 
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12.  EVALUATION  IN  A  BRANCH-AND-BOUND  METHOD 


The  other  element  needed  In  the  algorithm  is  an  efficient  method  of  bounding 
the  optimal  value  of  the  functional,  F  ,  from  below  In  any  partially  arbitrated 


problem.  Thus,  for  any  arbitration  a  in  the  arborescence  of  Figure  14 

* 

(a  -  0,1,  ....  A-Aq)  ,  we  know  that  there  will  be  a  best  value  for  F(a)  ■  F  (a) 


in  the  reduced  problem  governed  by  (22)  for  Che  already  arbitrated  values  of 
What  we  seek  is  an  evaluation,  E(a)  ,  of  this  optimum  satisfying 


(23)  E£al  <_  F*(a)  <_  E(a)  , 

one  which  can  be  found  without  calculating  all  of  complete  solutions  (terminal 

nodes)  which  are  descendants  of  node  a  .  (It  also  turns  out  to  be  easy  and 

worthwhile  to  find  an  upper  bound,  E(a)  ,  at  the  same  time,  as  indicated.)  These 

bounds  are  then  used  to  guide  the  exploration  of  the  arborescence  in  an  efficient 

A-Ao 

manner,  thus  avoiding  an  exhaustive  check  of  all  possible  2  solutions  to  the 

primitive  problem  (20). 

For  example,  at  the  parent  node  a  »  0  of  Figure  14,  we  have  the  original 
movable  activities  problem.  Call  F^  the  optimal  value  of  F  for  the 
corresponding  divisible  activities  problem  (1),  and  F°  the  critical  path  solution 
of  (1)  with  all  t^  *  0  ;  it  follows  that 


A  *  0 

t  <  F  (0)  <  + 


are  good  evaluators  in  the  sense  that  there  are  networks  for  which  one  (or  both) 
of  the  evaluators  will  be  attained  by  F*  . 

Similarly,  define  F^(a)  and  F^(a)  as  the  optimal  values  of  the  divisible 
activities  and  ordinary  critical-path  functionals,  respectively,  of  the  problem  of 
arbitration,  reduced  by  applying  (22).  Then 
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(25)  Eisi  -  F#(a)  F* (a)  <  F°(a)  +  J  T,  -  E(a) 

unarbitrated 
groups  d 

are  also  good  evaluators,  in  the  same  sense.  (25)  also  has  the  desirable  property 

that  at  a  terminal  node  of  Figure  14  (all  arcs  arbitrated),  the  bounds  coincide, 

* 

and  E  *  F  »  E  . 

Unfortunately,  using  the  divisible-activities  algorithm  is  very  inefficient. 
Not  only  may  several  simplex  iterations  be  necessary  to  solve  for  a  single  F^(a) 

but  columns  obtained  in  this  solution  cannot  be  used  for  another  evaluation 

# 

F  (a')  ,  since  they  may  not  be  feasible  in  the  new  subproblem  of  different 
arbitration. 

However,  we  note  that  repeated  computation  of  the  ordinary  critical  path 
solution  can  be  easily  done  for  any  arbitration,  finding  not  only  F°  ,  but  both 
the  early  jv^j-  an^  latc  jv^|  event  times  for  all  nodes.  This  observation  leads 
to  a  weaker,  but  much  rapid  determination  of  a  bound,  E(a)  . 
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13.  AN  EFFICIENT  LOWER  BOUND 

As  in  the  algorithm  of  Section  5,  define  the  slack  of  ij  for  any  ordinary 
critical  path  solution  to  be 


(26) 


s  .  _  VL  -  T 

ij  Vi  Vj  Xij 


Now  suppose  that  there  was  only  one  group,  ,  of  arcs  left  to  arbitrate. 


Def ining 


(27) 


S,  -  max  (s  )  , 
“  1J 


we  see  that  an  exact  solution  of  the  one  group  problem  is 


(28) 


E  =  F  +  max  (0,T1  -  S^)  =  F  ; 


in  other  words,  the  remaining  task  is  allocated  to  the  largest  slack  location  (or 
possibly  to  one  of  several  locations,  if  >  T^)  . 

When  there  is  more  than  one  group  to  be  arbitrated,  the  best  possible 
situation  is  to  pick  that  group  first  which  gives  the  largest  increase  in  ; 

this  increase  must  be  absorbed  in  any  case,  and  it  might  then  be  possible  to  tuck 
in  the  remaining  tasks  in  parallel  gaps. 

In  other  words,  at  arbitration  a  ,  solve  for  the  two  critical  path  solutions 
|v^(a)j-  ,  |v^(a)j  and  the  associated  duration  F^(a)  =*  v|j(a)  -  v^(a)  ■ 

v^(a)  -  v^(a)  ,  using  the  appropriate  .  Then 

(29)  E(a)  «  F^(a)  +  max  ("0  ,  max  (T^  -  Sd(a))"| 

I  unarbitrated  I 

{_  groups  d  J 


is  an  appropriate  evaluator. 
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Note  that  if  S,(a)  >  T,  for  all  remaining  groups,  we  cannot  necessarily 
d  —  d  ' 

conclude  that  the  problem  is  over,  since  the  largest  slacks  for  each  group  may  not 
be  simultaneously  available. 
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14.  IMPROVING  THE  UPPER  BOUND 

The  upper  evaluator,  E(a)  ,  given  in  (25)  can  easily  be  improved  by,  say, 
completing  the  arbitration  a  in  any  arbitrary  manner,  and  taking  the  resulting 
terminal  . 

A  simpler  method,  which  does  not  require  a  third  critical  path  solution,  is  to 
use  the  found  in  (27).  Suppose  we  pick  the  group  to  allocate  next  which  gives 

the  smallest  increase  -  S^  .  ^hen,  the  worst  case,  all  of  remaining  slacks 
might  have  decreased  to  value  zero,  and  all  remaining  unarbitrated  jobs  would  have 
to  be  added  in  series  to  the  critical  path.  Since  we  can  pick  this  first  group 
judiciously 

(30)  E(a)  -  F°(a)  +  £  -  max  min  [S,(a),T.]  . 

unarbitrated  unarbitrated  C 

groups  d  groups  d 

One  can  improve  this  bound  by  looking  ahead  several  steps,  but  it  is  probably  not 
worth  the  extra  effort. 
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15.  SELECTING  THE  ARBITRATION 

As  is  usual  in  branch-and-bound  algorithms  where  the  bounding  procedure  does 
not  lose  efficiency  when  skipping  among  the  different  candidates  in  the 
arborescence ,  we  recommend  that  the  arbitration  with  smallest  E(a)  be  separated 
next  during  the  algorithm. 

If,  at  any  stage,  we  discover  E(a^)  <  E^)  for  two  candidate  arbitrations, 
we  can  truncate  the  arborescence  at  node  ,  i.e.,  discontinue  further 
exploration  in  descendant  nodes.  This  illustrates  possible  value  of  carrying  upper 
bounds. 

Selecting  the  next  movable  task  to  be  arbitrated  should  be  done  adaptively, 
and  not  in  a  fixed,  predetermined  order.  Limited  computational  experience  suggests 
that  the  following  rule  is  quite  efficient: 

Select  as  the  next  movable  task  to  be  separated  from 

arbitration  a  that  unarbitrated  task  d  which  maximizes 

(31) 

T,  -  S,  (a)  .  If  there  is  a  tie  (or  if  all  T,  -  S  (a)  are 
ad  a  a 

zero),  select  the  task  d  with  largest  . 


In  other  words,  it  seems  efficient  to  decide  early  where  to  put  the  jobs  with  the 
most  excess  duration;  or,  if  there  is  a  tie,  to  decide  the  job  with  largest 
duration . 


Furthermore,  we  note  when  the  separation  occurs,  we  actually  know  in  advance 

what  the  values  of  F^(a  +  1)  will  be,  from  the  value  of  T,  of  the  {s..} 

d  ij 

belonging  to  this  group.  It  is  then  efficient  to  carry  out  the  evaluation  of  this 


separation  in  order  of  decreasing 


Often  the  upper  bound  from  one 


placement  of  in  a  large  slack  location  will  truncate  the  further  exploration 


of  the  arborescence  from  another  location  with  small  slack,  thus  saving  some 
computation.  This  observation  is  particularly  true  towards  the  bottom  of  the 
arborescen:e,  when  the  spread  between  upper  and  lowe *  evaluators  is  small. 
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16.  THE  MOVABLE  ACTIVITIES  ALGORITHM 

Figure  15  shows  the  flow  chart  of  the  final  algorithm.  Several  explanatory 
comments  are  perhaps  In  order. 

At  each  stage  of  the  algorithm,  a  candidate  list,  j?  ,  of  arbitrations  (nodes) 
is  kept;  these  are  merely  partially  arbitrated  solutions  which  have  not  been 
further  explored.  If  one  candidate  is  completely  dominated  by  another,  then  it  is 
eliminated  from  £  ,  and  the  arborescence  is  truncated,  as  described  above. 

Since  the  next  arbitration  in  £  to  be  elaborated  further  is  the  one  of 
current  smallest  £  ,  the  procedure  may  jump  between  nodes  of  different  order,  and 
the  resulting  sequence  will  not  follow  the  prior  numbering  of  the  nodes.  This  rule 
requires  more  computer  storage  than  other  methods,  such  as  "backtracking,"  but  it 
does  enable  one  to  deduce  that  an  optimal  solution  has  been  found,  whenever  an 
arbitration  of  maximal  order  is  reached  [1],  [8]  (assuming  all  ties  in  selecting 
JE  have  been  resolved)  . 

Various  other  heuristics,  perhaps  utilizing  the  upper  evaluator  and/or  the 
level  of  arbitration,  are  of  course  possible. 
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FIGURE  15: 


FLOW  DIAGRAM  OF  MOVABLE  ACTIVITIES  ALGORITHM 
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17 .  EXAMPLE 


Consider  Figure  16,  which  shows  the 


‘ij 


for  each  arc. 


four  movable  groups: 


Suppose  there  are 


A 

A 

A 

A 


1 

2 

3 

4 


{(1,2), (1,3), (1,4), (5 ,8), (6, 8), (7, 8)} 
[(2, 3), (3, 4), (5, 6), (7, 6)} 

{(2, 5), (3, 6), (4, 7)} 

{(3, 5), (4, 6)} 


5 

6 

4 

8 

23 


From  the  ordinary  critical  path  solution.  Figure  17,  we  find  that  the  slacks  are 
(sane  order  as  above) . 


{s.j  |  ij  e  Aj)  =  (0,0, 2, 4, 0,1}  1 

{s^  |  ij  e  A2)  =  {0,1, 2, 4}  2 

{ s±j  |  ij  e  A3}  -  {2,0,1}  2 

{s  .  |  ij  e  A  }  =  {5 , 3 }  3 

ij  4 


Ihu>  lot  the  patent,  node,  a  =  0  ,  we  conclude 

26  +  3  =  29  <  T*  t  26  +  23  -  5  =  4<*  . 


Since  is  largest  for  group  ,  we  arbitrate  this  first,  a  =  1  selects 

(3,5)  ,  a  =  2  selects  (4,6)  .  Note  that  we  have  immediately  the  't  l 
estimates  E( 1)  =  29  ,  E(2)  =  31  .  Following  cur  previous  rule,  wc  evaluate 
node  2  first,  cVt  lining 


E(2)  »  31  ,  E(2)  »  2 9  +  13  -  6  =  38 
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FIGURE  16:  MOVABLE  ACTIVITIES  EXAMPLE 


FIGURE  17:  ORDINARY  CRITICAL  PATH  SOLUTION 


FIGURE  18: 


OPTIMAL  SOLUTION' 
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from  Che  slacks  {0,3, 1,9, 0,6}  ,  {0,0, 7, 7}  ,  and  {7,5,6}  of  groups  numbers  one, 
two,  and  three  respectively. 

For  node  1  ,  we  obtain  slacks  {0,3, 5, 2,0, 4}  ,  {0,4, 0,5}  ,  and  {5,3,4}  for 
the  same  slacks,  giving  evaluations 

E(l)  =  29  +  1  -  30  ,  E(l)  =  29  +  15  -  5  -  39  . 

Since  node  1  has  the  lowest  E.  ,  we  arbitrate  that  node  on  group  2  ,  for 

which  T^  -  S^(l)  is  largest,  equal  to  1  .  This  creates  new  nodes  3,4,5, 

6  corresponding  to  selecting  (2,3)  ,  (3,4)  ,  (5,6)  ,  or  (7,6)  ,  respectively. 

We  can  get  estimates  of  =  35  ,  31  ,  35  ,  30  for  each  of  these  nodes  from  the 

last  step,  and  so  we  carry  out  the  evaluations  in  the  order  a  =  6  ,  4  ,  3  ,  5  or 

a  =  6  ,  4  ,  5  ,  3  . 

For  a  *  6  ,  the  slacks  are  {0,3, 1,3, 0,5}  and  {6,4,0}  for  the  remaining 

groups  numbers  one  and  three,  so  that 

E(6)  *=  30  and  E(6)  =  30  +  9  -  5  =  34  . 

From  this  upper  bound ,  we  can  immediately  truncate  nodes  3  and  5  and  not 
evaluate  them. 

For  a  =  4  ,  the  slacks  are  {0, 3 ,7 ,4 , 1 ,0}  and  {7,5,0}  ,  giving 

E(4)  =  31  and  E(4)  -  31  +  9  -  5  -  35  . 

The  list  J?  of  the  algorithm  now  contains  nodes  2,4,  and  6  of  the 

arborescence .  We  arbitrate  node  6  as  the  most  promising,  choosing  group  number 

one  since  it  has  largest  T,  (no  terms  T,  -  S,  positive).  The  result  is  the 

a  ad 

creation  of  nodes  7  ,  8  ,  9  ,  10  ,  11  ,  12  corresponding  to  choosing  arcs 

(1,2)  ,  (1,3)  ,  (1,4)  ,  (5,8)  ,  (6,8)  ,  and  (7,8)  ,  respectively.  The  preliminary 

estimates  of  £  are  35  ,  32  ,  34  ,  32  ,  35  ,  and  30  ,  so  that  the  evaluations 
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are  to  be  made  in  the  order  a  ■  12  ,  8  ,  10  ,  9  ,  7  ,  11  . 

For  a  *  12  ,  the  slacks  in  the  reamining  group  number  three  are  {6,4,0}  , 
so  that  there  are  two  locations  (2,5)  ,  (3,6)  in  which  the  four  units  of  work 
could  go.  Hence, 

E(12)  -  E(12)  =  30  . 

And  we  are  done,  because  of  the  look-ahead  property  of  our  evaluators,  and  the 
fact  that  no  remaining  £(a)  ,  a  e  £  ,  has  value  30  . 

Thus,  after  six  (earliest  and  latest)  critical  path  solutions,  we  have  the 
optimal  arbitrations 


{(3, 5), (7, 6), (7, 8), (2, 5)}  or  { (3,5) , (7 ,6) , (7 ,8) , (3,6) }  . 


18.  EXTENSIONS 


One  obvious  remark  is  that,  if  a  movable  activities  algorithm  is  available, 
then  one  can  use  it  (with  integer  modules)  to  give  a  good  starting  solution  to  the 
completely  divisible  activities  problem.  However,  as  mentioned  previously,  it 
does  not  seem  efficient  to  use  the  decomposition  algorithm  as  an  evaluator  routine 
for  the  movable  activities  problem. 

Another  straightforward  extension  is  to  a  cost-time  (CPM)  divisible  model. 

The  usual  CPM  algorithm  traces  out  the  project  cost  versus  duration  by  means  of  a 
parametric  change  in  dual  flow.  This  method  is  too  laborious  when  divisible 
activities  are  added;  however,  since  decomposition  is  being  used,  one  could 
minimize  total  project  cost,  adding  a  constraint  of  fixed  duration  to  the  master 
problem  (or  vice-versa) . 

Other  desirable  extensions  would  be  to  include  costs  in  a  model  with  movable 
activities.  For  example,  the  movable  tasks  could  be  of  fixed  duration,  but  there 
would  be  cost-time  tradeoffs  available  on  ordinary  activities;  a  CPM  computation 
would  be  used  for  the  evaluators  on  total  project  cost,  and  certain  nodes  of  the 
arborescence  would  probably  be  truncated  due  to  infeasibility  with  respect  to 
project  duration.  In  another  direction,  one  could  assign  certain  costs  to  placing 
a  movable  activity  between  a  certain  pair  of  nodes,  or  between  a  dummy  pair;  the 
resulting  optznal  insertion  model  could  be  used  to  solve  cost-time  critical  path 
problems  with  piecewise-linear ,  but  discontinuous  ar.i  nonconvex  costs.  The 
essential  difficulty  here  is  to  construct  efficient  evaluator  functions. 
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